Design of a Lambda system for population transfer in superconducting nanocircuits 
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The implementation ol a Lambda scheme in superconducting artificial atoms could allow detec- 
tion of stimulated Raman adiabatic passage (STIRAP) and other quantum manipulations in the 
microwave regime. However symmetries which on one hand protect the system against decoherence, 
yield selection rules which may cancel coupling to the pump external drive. The tradeoff between 
efficient coupling and decoherence due to broad-band colored Noise (BBCN), which is often the main 
source of decoherence is addressed, in the class of nanodevices based on the Cooper pair box (CPB) 
design. We study transfer efficiency by STIRAP, showing that substantial efficiency is achieved for 
off-symmetric bias only in the charge-phase regime. We find a number of results uniquely due to 
non-Markovianity of BBCN, namely: (a) the efficiency for STIRAP depends essentially on noise 
channels in the trapped subspace; (b) low-frequency fluctuations can be analyzed and represented 
as fictitious correlated fluctuations of the detunings of the external drives; (c) a simple figure of 
merit for design and operating prescriptions allowing the observation of STIRAP is proposed. The 
emerging physical picture also applies to other classes of coherent nanodevices subject to BBCN. 

PACS numbers: 03.67.Lx,85.25.-j, 03.65.Yz 
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I. INTRODUCTION 

The rapid technological progress in quantum-state en- 
gineering in superconducting nanodevices demands for 
the implementation of new advanced techniques of quan- 
tum control. STIRAP— is a powerful method in quan- 
tum optics, which is still largely unexplored in the solid- 
state realm. Using AC driving fields in A configuration 
(see Fig. [TJa) a quantum M > 2-state system is trapped 
into a subspace spanned by the two longest lived states. 
Control in this trapping subspace can be achieved by adi- 
abatic time evolution induced by properly crafted pulses, 
allowing for instance to prepare a given target state^i^. 
Adiabatic passage used in STIRAP guarantees highly ef- 
ficient and selective population transfer in atomic and 
molecular systemsii^. 

In the last few years it has been proposed that mul- 
tilevel quantum coherent effects^ could be observed in 
superconducting nanodevices, for instance electromag- 
netically induced transparency (EIT)^ or selective pop- 
ulation transfer by STIRAP^^"^. This would be im- 
portant both from a fundamental point of view, since 
coherent dynamics in multilevel atoms clearly displays 
beautiful interference phenomena^, and for applications. 
These include the implementation of microwave quan- 
tum switchesi^, the manipulation of solid-state qubit 
circuitsii and the fascinating perspectives of coupling 
strongly such nanodevices to electromagnetioi^ii^ or 
nanomechanical quantized modest. Very recently few 
experiments have demonstrated features of multilevel co- 
herence in such devices, as the Autler-Townes (AT) split- 
tin g^^d^ , KIT—, preparation and measurement of three- 



state superpositionsii, dynamical AT controU^ and co- 
herent population trappingi^. 

In all the above experiments, except that of Ref. [Tsl . 
the multilevel system was driven in the Ladder configu- 
ration^. Indeed in order to implement a A configuration 
the device Hamiltonian should be strongly asymmetric, 
which may be achieved by a proper external biasingi^— li^, 
otherwise selection rules prevent to drive efficiently the 
pump transition. However the longest decoherence times 
in quantum bits are achieved by biasing the devices at 
(or near) parity symmetry points. Hence difficulties in 
implementing a A configuration in superconducting nan- 
odevices raise a fundamental design issue. In particu- 
lar, low-frequency noise which is known to determine the 
performances of systems operated as quantum bits^^ and 
which is minimized at symmetry points, is shown in this 
work to play a major role in degradation of efficiency in 
STIRAP. So far the effect of decoherence in multilevel su- 
perconducting artificial atoms has been addressed using 
Markovian master equations. In this work we address de- 
coherence effects due to a solid-state environment where 
a strong non-Mar kovian noise component is also present. 
From the exquisite sensitivity of coherence to operating 
conditions, and to design parameters of the device we 
determine the prescriptions for the demonstration of a A 
scheme in realistic superconducting nanocircuits. 

We tackle this problem by a quantitative analysis of 
a class of superconducting nanocircuits, namely those 
based on the Cooper pair bo»2S (CPB, see Fig.H]). This 
is an important case-study encompassing several differ- 
ent coherent nanodevices which have already successfully 
implemented quantum bits^i"— . The emerging physical 
picture is even more general holding true for nanodevices 
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FIG. 1. (color online) (a) Three-level system driven with 
AC fields in A configuration, (b) The counterintuitive se- 
quence: the Stokes field is switched on before the pump field 
(here ^IqT = 20, r = 0.6 T). (c) Instantaneous eigenvalues 
{eo{t),e±(t)}, for S = 0, Sp = -0.2 r2o and k = 1. (d) Popu- 
lation histories pii{t) = \{i\tp(t))f ior ideal STIRAP (3 = 0): 
the system prepared in |0) follows the Hamiltonian along the 
eo adiabatic path yielding complete population transfer to |1). 



suffering mainly from the presence of low- frequency noise. 

The main message of this work is twofold. First we find 
that observation of STIRAP should be possible with de- 
vices fabricable at present days, provided that operating 
conditions and suitable design optimize the conflicting re- 
quirements of efficient coupling between states with (ap- 
proximately) the same parity and protection from low- 
frequency noise. Second, despite of the complicated mul- 
tilevel structure and of the many parameters involved, we 
show that the efficiency for STIRAP depends essentially 
on noise channels involving the trapping subspace, and 
determine a simple figure of merit for design and oper- 
ating prescriptions of devices allowing the observation of 
STIRAP. 

The paper is organized as follows. We introduce STI- 
RAP in Sec. Inland describe population transfer via adi- 
abatic and nonadiabatic patterns. In Sec. IIIII we discuss 
the implementation in a CPB, and introduce the model 
for broad-band colored noise (BBCN), extending to a A 
system the approach introduced in Ref s . [27l and \2&L which 
quantitatively explains qubit decoherence due to BBCN 
in superconducting qubits22r^. In Sec. II VI we present re- 
sults on the effects of the BBCN, focusing on the charge- 
phase regime of CPB's. In Sec. |V| we extend the above 
considerations to other regimes of the CPB and deter- 
mine the figure of merit characterizing optimal design 
and operating conditions. In Sec. IVII we compare the 



effects of dephasing with long memory time with Marko- 
vian dephasing, showing that in the former case driv- 
ing more strongly the system would improve efficiency. 
Tlierefore STIRAP could in principle discriminate be- 
tween different dynamic characteristics of decoherence 
sources in superconducting nanocircuits. Conclusions are 
drawn in Sec. IVIII 



II. COHERENT POPULATION TRANSFER IN 
THREE-LEVEL ATOMS 

A. Dark state and STIRAP 

In quantum optics STIRAP is based on a A configura- 
tion (Fig. [1^) of two hyperfine ground states |0) and |1) 
and an excited state |2), with energies Eq = 0, Ei and 
E2 respectively. The system is operated by two classical 
laser fields^"— , the Stokes laser = fls cos ujst and the 
pump laser Qq2 = fip coswpt, each being nearly resonant 
with the corresponding transition. The effective Hamil- 
tonian is conveniently written in a doubly rotating frame 
at the angular frequencies Wfc, where k = p,s refer to 
pump and Stokes. In the rotating- wave approximation 
(RWA) the effective Hamiltonian isiiiiM 



i7 = <5|l)(l|+<5,|2)(2| + hf|2)(l| + ^|2)(0|+h.c. (1) 



Here 5s = E2 ^ Ei — and dp ~ E2 — Eq — ujp are 
the single-photon detunings and we introduced the two- 
photon detuning 5 = dp — 6s. Both the detunings and the 
Rabi frequencies flk can be functions of time. At two- 
photon resonance, (5 = 0, the Hamiltonian ([1]) has a zero 
energy instantaneous eigenvalue eo = (Fig. ^) whose 
eigenstate is a "dark state" , 



1^) 



»,|o) - r!p|i) 



(2) 



and two other eigenstates |±) with nonzero eigenvalues 

1 



flp whose form can be found 



analytically^. If the system is in the dark state, the 
population is trapped in the two lowest diabatic states 
{|0), |1)}. This is due to destructive interference of the 
two fields: despite excitation by the lasers the state |2) is 
never populated and no radiative decay can be detected. 

By slowly varying the coupling strengths, fts{t) and 
rip(t), the dark state can be rotated adiabatically in the 
subspace spanned by |1) and |0). In particular STIRAP 
yields complete coherent population transfer |0) — )• |1) as 
foUowsi: the system is prepared in |0), which coincides 
with the dark state for fig ~ Qp ~ 0; then fl^ is slowly 
switched on; after a delay r also Op is slowly switched 
on; at this stage is slowly switched off and the dark 
state now coincides with |1); finally the protocol ends by 
switching off f2p, achieving complete population transfer 
(Fig. [lid). Notice that in STIRAP population trans- 
fer is achieved by a "counterintuitive" pulse sequence 
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FIG. 2. (color online) Sensitivity to detunings of the efficiency 
of STIRAP. Here QqT = 20, r = T/2 and k. = 1. In the 
white zone the efficiency is larger than 90%. Efficiency is 
very sensitive to a nonzero two-photon detuning 5, and much 
less sensitive to (5p 7^ (notice the different scale of the axis). 
The two lines on the plot represent correlated stray detunings 
induced in the CPB by charge noise for two different values 
of qg = 0.47, 0.49 (see Sec HITDl and HV^ . 



(Fig. [T]b), which has several advantages'^. First the 
excited state |2), which may undergo strong spontaneous 
decay deteriorating the transfer efficiency, is never popu- 
lated during STIRAP. Moreover, provided adiabaticity is 
preserved, STIRAP is insensitive to many details of the 
protocol, as the precise timing of the operations, a prop- 
erty which makes it attractive for implementing fault- 
tolerant quantum gatesiii^. 



B. Sensitivity to parameters 

Adiabaticity is important to achieve high efficiency 
since non adiabatic effects trigger unwanted transitions 
detrapping the system from the dark state. A neces- 
sary condition for adiabaticity isi KOpfi^ — Up(ls)/{^p + 
^s)l ^ k± ^ £o| which suggests that large enough Rabi 
peak angular frequencies f2fe are needed, in order to de- 
termine a large AT splitting of the instantaneous eigen- 
states. We let Qpit) = flo f[{t - r)/T] and fi^W = 
kHq f[{t+T)/T] . A positive delay r implements the coun- 
terintuitive sequence. For Gaussian pulses, f{x) = e^^ , 
the choice r > (V2 - 1)T and ^qT > 10 yields efH- 
cient population transfer—. As we discuss later, in su- 
perconducting nanocircuits the pump peak Rabi angular 
frequency cannot be very large, and T is limited by 
decoherencc. We found a good tradeoff for HqT = 15 and 
a delay r — 0.5 T, which turns out to be a satisfactory 
choice when fluctuations of parameters are considered. 
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FIG. 3. (color online) Non-ideal STIRAP {QoT = 20, r = 
0.6 T) with S ^ 0, shows different classes of patterns of in- 
stantaneous eigenstates. (a) Top: instantaneous eigenstates 



for 6 = 0.25 iln and Sr, 



-0.4 f2o. These LZ patterns with a 



single avoided crossing during the pump induced EIT phase 
result from the effect of low-frequency charge noise in CPB's 
(or flux noise in flux qubits). Bottom: population histories of 
the diabatic states, (b) Top: generic LZ pattern {5 = 0.5 fio, 
5p — 0.5 fio); bottom: population histories. 



and which we use unless otherwise specified. 

Non-zero detunings Ss and 6p modify the whole adi- 
abatic picture of STIRAP and may strongly affect the 
transfer efficiency. The crucial parameter is the two- 
photon detuning since for (5 7^ the dark state ^ is 
not anymore an instantaneous eigenstate and there is 
no adiabatic connection from the initial to the target 
state. As a consequence the efficiency is very sensitive 
to fluctuations of the two-photon detuning (5, whereas 
large single-photon detunin gs 5 r, are tolerable (see Fig. [2] 
and the discussion in Ref. l34l ). For 6=^0 the simple 
picture of adiabatic passage is not valid anymore and 
qualitatively new phenomena occur enriching the phys- 
ical scenario. In particular non-ideal STIRAP may still 
take place via non-adiabatic transitions between adia- 
batic states. For small values of S, narrow avoided cross- 
ings between the instantaneous eigenvalues occur and the 
population is transferred by Landau-Zener (LZ) tunnel- 
ing^ (see Fig. [3|) . For increasing S the transfer efficiency 
is reduced and in general the excited state |2) is popu- 
lated during the protocol. 

It is worth stressing the importance of correlations 
between detunings. Indeed it is well known in atomic 
physics^ that if Ss and Sp are correlated as to nearly 
preserve two-photon resonance, still a large transfer effi- 
ciency is obtained. In superconducting nanodevices cor- 
relations of other nature may arise between effective fluc- 
tuations of S and Sp, induced by solid-state noise. These 
correlations are represented by the lines in Fig. [51 and 
determine the typical pattern for population transfer via 
LZ processes^ of Fig. [3^. Notice that in this case rela- 
tively large single-photon detunings Sp ^ 25 S still allow 
coherent population transfer. This has important conse- 
quences in coherent nanodevices where fluctuations may 
produce large detunings Sp. 
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FIG. 4. In the CPB design the state of the superconducting 
island is a superpositions of states with well defined number 
n of extra Cooper pairs. The device is biased by the gate 
voltage Vg determining the operating point of qg = Vg/2eCg; 
control is operated by an ac component of Vg. Charge fluctua- 
tions are equivalent to voltage fluctuations 5Vx- The effective 
Josephson energy can be tuned via the flux "I>g of the magnetic 
field threading the loop, Ej — Ej{^g). 



III. STIRAP IN THE COOPER PAIR BOX 

A. Implementation of the A system 

The CPB^ is a superconducting loop interrupted by 
two adjacent small Josephson junctions (energy Ej/2) 
defining a superconducting island (Fig. |4|). The total 
capacitance C gives the charging energy Ec ~ (2e)^/2C. 
The electrostatic energy is modulated by a gate voltage 
Vg, connected to the island via a capacitance Cg C. 
The Hamiltonian reads 



Ec{n-qg)'\n){n\~^{\n){n+l\+h.c.) (3) 



where {\n), n €] — oo,oo[} are cigenstatcs of the num- 
ber operator h of extra Cooper pairs in the island. We 
have defined the reduced gate charge Qg = CgVg/{2e) po- 
larizing the island. The spectrum can be modified by 
choosing a specific bias qg (Fig. [5]). 

The parametric dependence of on qg defines a port 
allowing for external control of the system: by adding 
an ac microwave component qg ^ qg + qdt), shaped in 
suitable pulses, arbitrary rotations of the quantum state 
have been demonstrated^. 

In the basis of the eigenvectors {\(j)i{qg)), i — 0,1,2} 
of H{){qg) the driven Hamiltonian reads 

H{f) = Y,EM^){^^\ + A{t) Y^n,, (0,1 (4) 

i ij 

where = (0i|n|0j) and the control field is A{t) = 
-2Ecqc{t). For STIRAP we let A{t) = As{t)cosLOst + 
Ap{t) cos uipt. We then transform the Hamiltonian to the 
doubly rotating frame, and retain only slowly varying 
terms, which yields the RWA (see App. lA By project- 
ing onto the three lowest levels, i,j = 0, 1,2, we finally 
obtain an effective Hamiltonian H implementing the A 
configuration of Eq. Q , with the definitions 
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(5) 



Therefore, h enters the peak Rabi angular frequencies, as 
the electric dipole does in atoms. The CPB is a tunable 
atom since the parametric dependence on qg (see Eq.®) 
affects "diabatic states", eigenenergies, and matrix ele- 
ments of h (see Fig. [5]). Therefore, detunings and peak 
Rabi frequencies in H Eq.([T]) depend on qg. 

Several superconducting qubits are based on the CPB. 
From the point of view of the model they differ for the 
values of the parameter J = Ej/Eq. Computational 
states, which are eigenstates of Hq, are superpositions 
of a number of "charge states" |n), increasing with J, 
and therefore these devices have very different energy 
spectra. Coherent dynamics has been observed in the 
charge regime^ii^ J ^ 1 , in the charge-phase regime^^i^ 
J ^ 1, and in the phase regime^^ J ^ 1 (from several 
tens up to several hundreds). Physically these devices 
greatly differ both in the design (size, on-chip readout 
scheme) and in characteristics (ease of coupling to control 
fields, resilience to noise), these features being crucial for 
functionality. Therefore the CPB allows for a thorough 
discussion of requirements to observe STIRAP in a wide 
class of nanodevices. 

Notice finally that, besides transitions to higher en- 
ergy levels, the external field coupling to artificial atoms 
may in principle trigger the — t- 1 transition and for 
qg 7^ 1/2 also provides a time-dependent diagonal con- 
tribution to the effective Hamiltonian. We have shown 
in previous works that this has no effect for CPB's in 
the charge-phase regime, due to the large anharmonic- 
ity of the spectrum, even in the presence of markovian 
noise^i^i^. Therefore we will safely study the effect of 
noise in the lowest three-level subspace. Leakage from 
this subspace is expected for J ':$> 1, but in this regime 
STIRAP does not occur even in the three-level approxi- 
mation, as we show later. 



B. Symmetries, decoherence, selection rules 

Tunability with qg has been exploited to find optimal 
points where qubit operations are well protected from 
low-frequency noise22i21. For instance for qg = 1/2 the 
Hamiltonian ([3]) is symmetric for charge-parity transfor- 
mations (App. I A 2p . Due to this fact CPB-based qubits 
biased at symmetry are well protected against external 
noise. This has allowed to obtain experimental dephas- 
ing times of several hundreds of nanoseconds in charge- 
phase devices^^i^ and ranging from T2* > 2 ^s J ~ 50 
in the phase regime2^ up to T2 ^ 0.1ms recently re- 
ported^. At the same time symmetry enforces a se- 
lection rule preventing transitions between states with 
the same charge-parity. In particular no2 vanishes at 
qg = 1/2 (see Fig. [SJd) therefore it is not possible to im- 
plement the A configuration of Eq.([T|), since fip = 0. In 
Refs.l^ and|^ it has been proposed to overcome this prob- 
lem by working slightly off-symmetry (sec Fig.[S]D,c), and 
it has been shown that the full multilevel structure of a 
CPB with Ej = Ec allows for coherent population trans- 
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obtained on a phenoinenological level by "quantizing" 
noise. This is obtained by letting SH = X h + Hji, where 
X is an environment operator and Hn describes the en- 
vironment alone, and suitable counterterma^. Marko- 
vian noise can then be studied by deriving a weak cou- 
pling quantum optical master equation (ME). However, 
noise in the solid state has large low-frequency compo- 
nents invalidating the ME. A multistage approach has 
been proposed^! where high and low-frequency noise are 
separated, the latter being approximated by a classical 
random field. Formally X — )■ Xf — 2Ec x{t) where Xf 
describes fast environmental quantum degrees of freedom 
and x{t) is a slow classical stochastic process. If we let 
Qx{t) = Qg + x{t) the Hamiltonian is written as 



FIG. 5. (color online) (a) Energy spectrum Ei of a charge- 
phase CPB for J — 1.32 (corresponds to the Quantroniunj^^) , 
relative to the ground state Eo — 0, vs. the bias qg. (b) Ma- 
trix elements of h involved in the A scheme vs. qg for J = 1.32; 
the element no2 vanishes at the symmetry point qg = 1/2; 
(c) Matrix elements vs J = Ej/Ec for qg = 0.48; notice 
that no2 is much smaller than other elements (it vanishes at 
qg — 1/2) and it has nonmonotonous behavior for increasing 
J. 



fer for qg « 0.47, in the presence of Markovian noise. 

We stress that protection from noise and selection rules 
are related since they both stem from charge-parity sym- 
metry. Notice that increasing J enforces the (approxi- 
mate) selection rule in a larger and larger neighborhood 
of the symmetry point qg = 1/2, since it makes less ef- 
fective symmetry breaking terms (asymmetric charging 
energy). For instance Fig. [5t shows that no2 at off- 
symmetry [qg = 0.48) eventually decreases for increas- 
ing J, making impossible the implementation of the A 
scheme. 



C. Model for charge noise 

In principle each port of the device also allows injec- 
tion of noise and provides channels for decoherence. The 
control port associated to qg couples to charge noise and 
in this paper we focus on it, since it is the main source of 
low-frequency noise in the CPB for the regimes in which 
STIRAP could be observed. The structure of the cou- 
pling to noise can be obtained by allowing for fiuctua- 
tions of the gate charge in the Hamiltonian Eq. ([3]) . Their 
physical origin, besides voltage fluctuations of the circuit 
have been recognized as the effect of switching impuri- 
^ie j29^40,4i located in the oxides or in the substrate close 
to the device. We let qg ^ qg + x, where x describes 
stray electrical polarization of the island, and write the 
resulting Hamiltonian as H = Ho{qg) + Hnwit) + SH. 
Here Hjiwit) is the control Hamiltonian in the RWA, 
Eq. ()A3|) . whereas 6H = —2Ec xn describes fiuctuations. 
The structure of coupling to a quantum environment is 



H = Ho[qg + x{t)] + HRw{t) +Xh + He 



(6) 



In many cases low- frequency noise has 1// spectrum 
and the leading contribution of the slow dynamics of 
x(t) is captured by a static-path approximation (SPA) 
i.e. approximating the stochastic process by a suitably 
distributed random variabl o^^'^° x. In this simpler sce- 
nario one should first calculate the reduced density ma- 
trix p(t\x) for a given stray bias x obtained by tracing out 
high-frequency (quantum) noise, and then average over 
the distribution p{x). In particular population histories 
are given by Pi{t) = J dxp{x) pii{t\x). Notice that for 
each realization x of the random variable the system is 
prepared and measured in the eigenbasis of Ho(qg + x), 
which is then conveniently used to represent pu. 

In the case of many weakly coupled noise sources p{x) 
is a Gaussian with standard deviation ax- The low- 
frequency noise affects the dynamics via fluctuations of 
energy it induces. This point of view provides a simple 
argument explaining why the symmetry point qg — 1/2 
is well protected. Indeed, since at this working point 
the energy splitting Ei depends only quadratically on 
the fiuctuations x, energy fluctuations are suppressed. 
Therefore superpositions of the two lowest energy levels 
keep coherence for a longer dephasing time, with only 
a power law suppression of the signal^L^. This case 
is referred as "quadratic noise" regime, to make a dis- 
tinction with "linear noise" conditions, occurring for off- 
symmetry bias, where energy fluctuations are linear in 
X yielding much stronger decoherence (Gaussian decay 
law). 

This approach has quantitatively explained the power 
law decoherence observed not only in CPB's^, but also 
in flux qubits^ and has allowed to find optimal operating 
point in ultrafast driven phase qubits^^. Recently it has 
been used to discuss properties of multiqubit systems^. 
The present extension to a A system of the approach of 
Ref. 123 enlightens the role of correlations between detun- 
ings, and provides a tool for optimal device design. 
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D. Effective model for low- frequency noise in A 
configuration 

In order to study STIRAP the Hamiltonian (|6]) is pro- 
jected onto the subspace spanned by the three lowest 
energy adiabatic eigenvectors of Ho[qg + x{t)]. In do- 
ing so we assume the adiabaticity of the dynamics in- 
duced by x{t), which aUows to neglect effects of the timc- 
dcpendcncc of the eigenvectors. Of course in the SPA 
adiabaticity of noise is automatically verified. The sys- 
tem plus drive Hamiltonian HQ[qg+x{t)]+Hfiw(t) in the 
rotated frame has the same structure of Eq. ([!]), but de- 
pends on the realization of the stochastic process. Fluc- 
tuations of the eigenenergies translate into fluctuations 
of the detunings (we let Eq = 0). In the SPA we have 

S{x) = Ei{qg + x)~ujp+u)s ; Sp{x) = E2{qg+x)-ujp (7) 

It is worth stressing that also the effective drive fluctu- 
ates, via the charge matrix clement, for instance flp = 
"■02 [^g + x{t)] Ap. Thus the effect of low- frequency noise 
in solid-state devices is conveniently recast in terms of 
sensitivity of the protocol to flctitious imperfections (in 
both phase and amplitude) of the drive. This allows to 
apply to solid state devices several results from the quan- 
tum optics realm. For instance, the known critical sensi- 
tivity to two-photon detuning, translates in the fact that 
the main flgurc to be minimized, in order to achieve effi- 
cient population transfer in nanodevices, are fluctuations 
of the lowest energy splitting. This is a quantity which 
is well characterized from the qubit dynamics^Sri^. 



IV. EFFECT OF BROADBAND COLORED 
NOISE ON STIRAP IN A CHARGE-PHASE CPB 

We now apply the above approach to analyze STIRAP 
in a CPB in the charge-phase regime Ej ~ Ec- An im- 
portant point is that while dephasing is minimized by 
operating at the symmetry point qg = 1/2, the selec- 
tion rule no2 ~ apparently prevents to implement STI- 
RAP. Therefore it has been proposed to operate slightly 
off-symmetryii^, where on the other hand decoherence 
due to low-frequency noise increases^. This opens the 
question of the tradeoff between efficient coupling of the 
driving fields and dephasing due to slow excitations in 
the solid-state. 

Since it is convenient to work with the largest possible 
pump Rabi peak frequency flo, wc will consider its value 
as a scale. For a given peak value of Ap{t) it can be 
estimated as ^Iq = il^j no2('j'g)/fioi(l/2) where flu is the 
maximal angular frequency for Rabi oscillations between 
the lowest doublet at the symmetry point, which is well 
characterized in experiments. We will use frequencies 
corresponding to vj^ ~ 600 MHz, which are in principle 
achievable^ even if there may be technical problems in 
specific devices. For the Quantronium at qg — 0.48 this 
would correspond to a maximum Vp = 55 MHz. 
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FIG. 6. (color online) Population histories in the Quantron- 
ium at qg = 0.475, averaged over fluctuations with ax ~ 0.004. 
Charge fluctuations determine anticorrelated stray detunings, 
Sp — —23.5 5. Drives are symmetrized, k = 1. The total time 
for the Quantronium corresponds to ~ 250 ns. The resulting 
efficiency is Pi{tf) = 0.77. 

Close to the symmetry point coupling of the field with 
the Stokes transition is larger. Therefore, wc could easily 
choose K = Vsjvp K, nY2l'nQ2 ^ 1- However, using larger 
values of docs not improve the transfer efficiency in 
CPB's^id^, therefore we will let k = 1 hereafter. 

A. Effects of low-frequency noise 

Low-frequency fluctuations x of the gate charge deter- 
mine non-exponential dephasing in qubits^i. They have 
been well characterized in the Quantronium by Ramsey 
interfcrometry at different bias points qg G [0.4, 0.5)^. 
Gate charge fluctuations (Jx are obtained by the mea- 
sured charging energy fluctuations oe = '^Eqctx ^ 
0.01 £'i(l/2); a flgurc which is independent on the bias, 
corresponding to = <je/{'^Ec) ~ 6-10~"^. These quan- 
tities are related to the integrated spectral density of the 
environmental, and for 1//" noise they also depend on 
details of the protocol as the total measurement time. 
Even if this dependence is only logarithmic, one can take 
advantage from the fact that measuring the final popula- 
tion in STIRAP requires a lower statistics than Ramsey 
fringes. Therefore, for our purposes lower values of are 
well reasonable and hereafter we use = 0.004. This 
is a realistic figure not only for the Quantronium but 
for the whole class of CPB-based devices, since charge 
noise is ultimately determined by material issues which 
are constantly under investigation. 

We consider STIRAP for the optimal conditions of 
nominal single and two-photon resonance, 5 ~ 5p = Q. 
According to Eq.(I7|) fluctuations x determine a distribu- 
tion of stray detunings. For small we can approximate 

5{;x) Aix + ]^Bi x^ 

where Ai ~ (dEi/dqg) = Ecai{qg,J) and Bi = 
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(d'^Ei/dqg) = Ec bi{qg, J). It is worth stressing that 
arbitrary small fluctuations determine S{x) ^ 0, there- 
fore STIRAP may occur only via non-adiabatic patterns. 
In the same way also Sp{x) depends on the derivatives 
{d^E2/dql'). 

In Fig. [6] we plot the populations histories Pi (<) aver- 
aged over the fluctuations of x. These induce correlated 
fluctuations of both detunings and couplings riij . Device 
and bias parameters correspond to a Quantronium bi- 
ased slightly off-symmetry. It is shown that low-energy 
fluctuations determine a ~ 20% efficiency loss despite of 
the fact that protection from noise is greatly reduced. 
This is an interesting figure for superconducting nanodc- 
vices if we compare with the observed coherent popula- 
tion trapping of ~ 60% recently measured in phase-type 
devices Moreover the population of the intermediate 
level is very small during the whole procedure, fulfilling 
the requirements for coherent population transfer. 

Such numerical evaluations are performed by using a 
4-th order integration Runge-Kutta method for the solu- 
tion of the ordinary differential equations. Convergence 
was tested down to a relative error lower than 10"'^, by 
adjusting both the integration step and the number of 
event series. Integration over fluctuations was performed 
by a Montecarlo approach with up to 5000 samples in 
order to attain a relative error smaller than 10"'^. 

Notice that the detunings depend on a single random 
variable x, therefore their fluctuations are correlated. In 
particular, charge noise determines anticorrelated fluctu- 
ations of effective detunings in CPB's, as it is clear from 
the spectrum (Fig.[S]a). This implies that non- ideal STI- 
RAP may occur only via the typical LZ patterns^ shown 
in Fig. |31a. 

Notice that in the regime of Fig. [5] fluctuations of the 
couplings riij could have been neglected. Indeed, they 
can be estimated from Fig. [5]b. For instance, for J ^ 1 
and Qg < 0.49 fluctuations of the amplitude of the pump 
pulse are Up ~ ao2 0'x^o, where 092 = dnQ2/dqg, there- 
fore (Tp <C r^o- Numerical results (Figs. ISIlOp actually 
confirm that fluctuations riij yield at most corrections, 
and moreover when they are appreciable STIRAP does 
not work due to the combined effect of high-frequency 
noise (see ^IVCp . 

The above observation implies that for practical pur- 
poses efficiency can be discussed entirely in terms of 
the sensitivity to detunings^. Diagrams in the {5, dp) 
plane (Fig. [2]) can be used to understand the effect of 
low-frequency noise. Correlated stray detunings in the 
CPB are there represented by lines in the 6 ~ dp plane, 
which are straight lines for linear noise. We draw the seg- 
ment corresponding to fluctuations x € [—C7x,c7x] (here 
ax = 0.004), for each bias point {qg = 0.47,0.49 are 
shown, the slope increasing by approaching qg = 1/2). If 
segments lie inside the light zone the efficiency is large. It 
is seen that efficient STIRAP requires small ffuctuations 
\S\ < 0.1 fig but the large anticorrelated \dp\ < 2.5 fto is 
tolerable. 
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FIG. 7. (color online) Population histories paify in the pres- 
ence of high-frequency noise (Pi is in absence), at resonance 
(5 = 5p = 0, left panel) and for finite anticorrelated detunings 
(5 = 0.05, 5p = —255, right panel), for k = \. Solid lines are 
obtained by inserting only 701 = l/Ti in Eq.([8]), describing 
relaxation 1 ^> only, and the associated secular dephas- 
ing. We have chosen a rather large 701 /fio = 0.01 to empha- 
size the effect. Dashed lines take into account all the other 
low-temperature emission and drive-induced absorption chan- 
nels (the chosen rates overestimate these processes), which are 
seen to have a limited impact on the efficiency. Physical scales 
for fio = 3.46 X 10*rad/s (the value we use for the Quantro- 
nium at qg = 0.48) are T ~ 43 ns, for the overall protocol 
Tt f» 290 ns and the chosen Ti^Tt. 

B. Effect of high-frequency noise 

High-frequency noise is studied by solving the 
quantum-optical ME in the rotating framed p = 
j^[p,H] — Dp, where p is the density matrix and H is 
the Hamiltonian ([T|). The structure of the dissipator Dp 
in the basis of the diabatic states {|0i)} reads^ 

(Dp)ij = ^ Ptj - Srj Pkkltk + (1 - 5tj)ltjPt3 (8) 

The first two terms describe emission and absorption of 
energy and the associated secular dephasing: 7^^ = ^j^i 
are transition rates between diabatic states, and 7^ = 
J2k^i^ki are the total decay rates of states \(j)i). At low 
temperature in an undriven system only rates of sponta- 
neous emission between diabatic states are non negligi- 
ble. In AC driven systems rates describing environment- 
assisted absorption are also nonzero, when the corre- 
sponding field is switched on^. Finally the dissipator 
may include pure dephasing rates jij = "fji. 

In quantum optical systems STIRAP connects two 
ground states, 701 = 710 = 0. Therefore as long as pop- 
ulation in |(/)2) is small all the transition rates act on de- 
populated states, and it is known that they practically do 
not affect population transfer. Instead in superconduct- 
ing nanocircuits the decay channel 701 is active. There- 
fore, we expect that 701 is the main source of efficiency 
loss due to processes involving energy exchange with the 
environment. This is indeed the qualitative conclusion 
suggested by the results in Fig. [T] 

To clarify the physical picture in Fig. [7] we study sep- 
arately the impact of adding decay channels. First we 
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consider only spontaneous deeay in the first doublet. We 
take 7oi/f2o ~ 0.01, which is a rather large value used 
to emphasize the effects and we study population histo- 
ries (solid lines We find pii{t) w Pi(t) e-T"i(*-*'\ 
where Pi{t) is the population in absence of noise, there- 
fore this channel mainly determines the simple popula- 
tion loss 1 — >■ when the target state is populated. It 
also determines a nonvanishing population P22 7^ which 
indicates detrapping from the dark state due to loss of 
coherence. 

Adding all the other decay channels (dashed lines) 
produces minor modifications of this picture (Fig. [7] 
left panel) for fields at resonance. No modification at 
all occurs for nonvanishing detunings, mimicking low- 
frequency fluctuations. The reason is that in this latter 
case 7oi does not determine substantial detrapping, and 
population of \(f>2)- 

In detail results of Fig. [7] were obtained by using 
rates for the other decay channels which overestimate 
unwanted processes, namely 712 = 2701, whereas 702 = 
0.2 701 (accounting for the suppression by selection rules). 
Notice that we do not take into account the fact that 
these emission rates become smaller when the drive 
amplitudes ftkit) are large enough, as resulting from 
the generalized (Bloch-Redfield) ME for AC driven sys- 
tems undergoing Rabi oscillations^. This latter ap- 
proach shows that also in the weak damping Rabi regime 
{Trtk{t) ^ 1 and 6k ^ ^k{t)) field- induced absorption 
sets in, even at low temperatures. We take into ac- 
count this channel phenomenologically, letting 721 (i) = 

712/4 [1 - 5,J{y^Sf+iW))?9mt)n where g(x) « 1 
only for x 3> 1 accounts for the requirement that field in- 
duced processes set in for underdamped Rabi oscillations. 
We used a similar expression for 720 (i)- 

Notice that while secular dephasing is taken into ac- 
count in Fig. [7] we did not include Markovian pure 
dephasing rates, jij = 0. Indeed, we argue that 
pure dephasing comes mainly from low-frequency (non- 
Markovian) noise accounted for by classical fiuctuations 
of X. In the next section we study the combined effect of 
high and low-frequency noise. We will discuss different 
models of pure dephasing in Sec. I VI I We finally mention 
that, for charge-phase CPB it has been shown that op- 
erating at Qg = 0.48 already provides sufficient coupling 
no2 to observe STIRAP in the presence of the Markovian 
component of noised. 



C. Combined effect of low and high-frequency noise 

The main conclusion of the last two sections is that 
the leading effects reducing coherent population transfer 
in nanodevices essentially involve decoherence of the first 
doublet. Another detrimental effect is that coupling to 
the pump pulse may be too weak due to (approximate) 
parity selection rules at (near) the symmetry point. 

With this in mind we investigate the interplay of low 
and high-frequency fluctuations, for qg < 1/2 in a charge- 
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FIG. 8. (color online) Efficiency of STIRAP Pi(f/) as a 
function of the bias qg in the presence of low-frequency and 
BBCN for the Quantronium [Ej/Ec = 1.32). Here ^oT = 15 
~ 0.004, 1/01 ~ 600MHz. Upper curves show effects of 
low-frequency noise, whereas the lower curve (black dashed) 
includes also high-frequency noise. Low- frequency noise is 
analyzed by adding different components, namely linear and 
quadratic correlated fluctuations of detunings (red solid curve, 
and red squares), linear and quadratic fluctuations of no2 
(blue solid curve, and blue dots). For off-symmetry bias 
{qg < 0.9), only linear detuning noise is important (see App.lBl 
for the behavior near qg < 1/2). In the inset the population 
P2(0) at intermediate times is shown. 



phase CPB. Indeed we will argue in the next section that 
STIRAP can be observed only in this regime. Here we 
consider a case-study device as the Quantronium, where 
noise in the first doublet has been well characterized. 
We take the value Ti = 1 /is which is achievable in the 
class of CPB devices at Qg = 1/2, and neglect its weak 
dependence on the bias^. 

Results are summarized in Fig. [8] where the efficiency 
is plotted against the bias qg, showing the impact of 
adding various low-frequency and high-frequency deco- 
herence channels. Curves refer to the same floT = 15, 
which guarantees adiabaticity for ideal STIRAP. It is ap- 
parent the different behavior sufficiently far {qg < 0.49) 
and close {qg ~ 0.5) to the symmetry point. 

For off-symmetry bias it is possible to observe STIRAP 
despite of the reduced protection from low-frequency 
noise. In this regime low-frequency noise is the main 
source of efficiency loss allowing a population transfer 
close to ~ 80%. Notice that only linear fiuctuations of 
the detunings are important: indeed Fig. |S] shows that 
accounting for the whole structure of low-energy fluctua- 
tions yields basically the same result, as the dependence 
of Eio and no2 on qg < 0.49 would a priori suggest. Effi- 
ciency is reduced to ~ 70% when also effects of high- 
frequency noise are taken into account (black dashed 
curve Fig. [S]). 
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Instead by approaching qg = 1/2, while low-frequency 
fluctuations would still allow for some population trans- 
fer, the interplay with high-frequency noise, mainly due 
spontaneous decay \(f>)i — ^ {(j^o), leads to the suppression 
of the efficiency (solid curve in Fig. |5]). Actually, in this 
regime the description of the effect of low-frequency fluc- 
tuations is more complicated spoiling the simple picture 
based on sensitivity to detunings. We discuss in App. |B] 
the whole information contained in Fig. [51 The main 
point is the observation that even if the device is biased 
at Qg = 1/2, fluctuations still allow for a nonvanishing 
pump coupling despite of the parity selection rule. The 
reason why high-frequency noise suppresses the efficiency 
is understood by recalling that STIRAP requires large 
pulse area, ^p{qg + x)T ^ 10. Since close enough to 
the symmetry point J7p becomes small, larger and larger 
T are needed which eventually exceed by far Ti. This 
mechanism explains the fact that the loss of efficiency 
due to high-frequency noise (Figs. I8ll0p appears to de- 
pend strongly on qg, even if we neglected the (in any case 
weak) dependence of the rates on the bias. 

The population P2 (0) of the intermediate level during 
the adiabatic passage phase remains small in the presence 
of BBCN (inset of Fig. [8]) . This is an essential require- 
ment for success and applications of the protocol, and 
completes the statement that STIRAP should be observ- 
able in charge-phase CPB's. 

Finally we mention that working with larger asymme- 
try, while providing a stronger pump coupling (see Fig.[S]) 
enhances the effect of low-frequency noise, reducing the 
overall efficiency. This is apparent from the trend in 
Fig. [8] (see also Fig. [9]), indicating that optimization of 
strong enough pump coupling and protection from low- 
frequency noise is a key issue for the implementation of 
a A system. 



V. OPTIMAL DESIGN OF THE DEVICE 

Efficiency of population transfer may be improved by 
optimizing the parameters of the protocol. In the last sec- 
tion we have shown that, due to the combined effect of 
the approximate symmetry and of spontaneous decay, ef- 
ficiency is large enough only if the device is biased slightly 
away from the symmetry point. In this section, we ar- 
gue that in this regime one should mainly optimize the 
tradeoff between coupling of the pump pulse and energy 
fluctuations of the lowest doublet of the device, due to 
low-frequency noise. Indeed the relevant figure of merit 
turns out to be 



2Ec{no2) ^ 
oc 



max 
P 



OS 
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(9) 



where os ^ y ^lO^ + ^B'^o'^ are the fluctuations of S{x). 
Efficient population transfer requires large enough values 
of this quantity, which depends on qg and on J, via (7102) 
and {Ai,Bi). This allows to choose convenient design 
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FIG. 9. (color online) The figure of merit Qp/as is plotted in 
the {qg, Ej/Ec) plane. We have chosen ~ 0.004 and f^o 
produced by an external field, which would determine Rabi 
oscillations with vr = 600 MHz in the first doublet. The 
analysis is valid far enough from the charge-parity symmetry 
point, which is not an interesting regime since the efficiency is 
suppressed. Dashed lines correspond to the values of Ej/Ec 
checked in this paper (Figs. 151 and 1101) . 



and operating conditions (see FiglH]). 

The criterion is clear from heuristic grounds, but can 
also be justified starting from an estimate of the linewidth 
for population transfer at finite 5. To this end we general- 
ize an argument given by Vitanov et ali^. They noticed 
that even if states of the adiabatic basis {\D), |±)} are 
not anymore instantaneous eigenstates, still \D) provides 
a connection between the diabatic states |0o) E^nd l^i). 
Then it is argued that efficiency loss depends on processes 
triggering transitions from \D) to |±). These are due to 
non vanishing off diagonal entries of the Hamiltonian in 
the adiabatic basis, which are proportional to 5. There- 
fore, if (5 ^ min|e±| population transfer does not oc- 
cur. This condition implies that for <5p = the linewidth 
scales linearly with the amplitude of the fields2i, 5i = 

^i''') \J i^p^"'^)'^ + (fi^°^)^- In our case stray detunings 
6{x) and 5p{x) ^ are anticorrelated, therefore leakage 
from \D) occurs during the pump phases (see Fig. Oa). 
Moreover, from Fig. [2] we see that it is substantial only 
when dp > 2D,o. In this regime the relevant condition 
Si = |e_| is an equation whose solution can be still 
written as 61 « (i'(T, k) f2™°^. Asking that fluctuations 
of 5 do not destroy the efficiency means that we need 
as ^ 5i. Therefore, we need large values of the parame- 
ter Si /(7s oc rt™°-^ /as, which justifies the figure of merit 
defined in Eq.(|9l). Our derivation does not take into ac- 
count fluctuations of the matrix elements, since they are 
negligible in the regime where STIRAP could work. For 
the same reason we did not include in Fig. [HI the region 
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cicncy {J — 0.7 in FiglTUk and J = 1 in FiglTUb). Instead 
for larger values of J, as in the Transmon^^ design, pump 
coupling is insufficient even if protection against noise 
is much better. In the opposite limit of charge qubits 
Ej /Ec ^ 1 , the efficiency is also small because of both 
small coupling and reduced protection from noise. This 
latter strongly suppresses population transfer also for a 
bias qg far-off symmetry, despite the coupling to the field 
increases. Notice that the dependence of the efficiency 
on the parameter J, besides providing prescriptions for 
the fabrication, can also be checked by on-chip tuning of 
Ej via an external magnetic flux (see Fig. 2]). 
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FIG. 10. (color online) Efficiency Piitj) vs. bias qg for J = 
0.7, 1,2 (see Fig. [9]). Parameters are the same as in Fig. [S] 
where Ti = 1000ns (black short-dashed curves). In the two 
upper panels efficiency for smaller Ti = 500 ns (gray short- 
dashed curve) is also shown. 



near qg = 1/2, since in this regime STIRAP is in any 
case prevented by spontaneous decay, due to the too low 
achievable values of fioT'i- 

We check the optimization suggested by FiglHlby look- 
ing at STIRAP for different values of J. It is seen that 
proper fabrication parameters allows to obtain larger effi- 



VI. COMPARISON OF DIFFERENT 
MECHANISMS OF DEPHASING 

Studying low-frequency noise in nanodevices by a non- 
Markovian model is necessary to explain quantitatively 
striking experimental features observed in quantum bits, 
as the peculiar non-exponential initial decoherence^"— . 
Moreover, this approach provides valuable additional in- 
formation as relations between effects of noise for differ- 
ent bias point^"— and different device design^, which 
are uniquely explained by the parametric dependence of 
the energy spectrum. We stress that such a picture is 
entirely due to non-Markovianity of BBCN. 

Therefore, this work complements previous studies in 
the quantum optics realm where typically the Marko- 
vian ME is used. In this latter approach pure dephasing 
is studied by considering only nonvanishing dephasing 
rates ^ij in the dissipator Eq.®, instead of the static 
fluctuations considered in this paper. Pure dephasing in 
the Markovian ME was studied by Ivanov et al42., who 
derived an adiabatic solution of the Liouvillc equation 
interpolating between the coherent and the incoherent 
limit. They predicted striking behaviors as a function of 
the control parameters, deriving several analytic results, 
which have been numerically checked. In particular for 
Gaussian pulses, populations at the end of the protocol 
were found to be 



1 



Pii(oo) = o + 



-37oiTV(8r) 



POO(OO) = P22(oo) = \- \^~ 



■37oiT^/(8r) 



(10) 



Notice that in this approximation the efficiency is de- 
termined by the dephasing rate of the lowest doublet 
only, 7oi. The conclusion that other dephasing channels 
are less relevant (actually for 712,702 ^ 701 some de- 
pendence appears in the numerical solutions of the ME) 
agrees qualitatively with our results with the static fluc- 
tuator model. On the contrary, the other striking feature 
of Eq.dTU]), namely that losses due to dephasing are in- 
dependent on the peak Rabi frequencies, does not hold 
for low-frequency noise. Following Ref. we plot in 
Fig. [TT] populations pii[oo) obtained numerically from 
the Markovian ME, using fixed T = = I/701, for in- 
creasing pulse amplitude fio- We compare them with the 
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FIG. 11. (color online) Efficiency of STIRAP (final popula- 
tions) as a function of the drive amplitudes fio. We compare 
the case of Markovian (pa) pure dephasing^ with the non- 
Markovian (Pi) model studied here. In both cases we let 
T = T2 = 57 ns, which for non-Markovian noise is obtained 
by taking ~ 0.004 in a device with J — 1.32 at qg — 0.48. 
It is seen that the effects of non-Markovian dephasing can be 
attenuated and suppressed by using larger Qo, whereas for 
Markovian noise STIRAP, when effective, does not depend 
on Qo- 



populations Pi (00) for the BBCN non Markovian model, 
where linear fluctuations of the detunings are considered 
such that ax — \/2l{A{T^), which yield the same in 
the qubit dynamics. It is seen that efficiency for BBCN 
depends on Oq ^nd improves for increasing values. 

Dependence on llo is a natural consequence of non- 
ideal STIRAP occurring via LZ patterns determined by 
low-frequency noise. Markovian noise cannot account for 
this scenario. The situation here is reminiscent of dy- 
namical decoupling^ which eliminates dephasing for 1// 
noise sources, as the effect of a strong continuous AC 
fields also does^. 

Another difference between Markovian and non- 
Markovian dephasing is that this latter practically does 
not populate the intermediate level |(/)2), although it de- 
creases the transfer efficiency. This is another indication 
of the reduced sensitivity of the protocol to low-frequency 
noise. On the contrary, sensitivity of P22 to Markovian 
noise is substantial and could give direct informations on 
7oi, as seen from Eq. lllOp . This observation is reminiscent 
of the proposal of Ref y of using EIT to probe decoherence 
of a phase qubit based on a SQUID nanodevice. Having 
in mind realistic noise spectra it is likely that the contri- 
bution of intermediate frequencies may determine effects 
similar to Markovian dephasing. Therefore, cross check- 
ing measurement of decoherence of two and three-level 
dynamics could give valuable spectral-resolved informa- 
tion on the environment. 

We stress the striking implication of non-Markovianity 



of the noise, namely correlations between fluctuations 
of the detunings, entirely determined by the paramet- 
ric dependence of the energy spectrum. Effects of time- 
correlated (Ornstein Uhlembeck) phase noise in opti- 
cal systems were studied by Monte Carlo simulations in 
Ref. l36l where the regime of partially correlated (5p and 
bs was addressed. In nanodeviccs the situation is differ- 
ent since we have strongly anticorrelated (or correlated) 
stray detunings 8 and 5p. It would be interesting to in- 
vestigate dynamic phase diffusion also in this case. 



VII. CONCLUSIONS 

In this paper we have studied the combined effect of 
low-frequency and high-frequency charge noise on the 
coherence of a CPB operated as a three-level artificial 
atom in Lambda configuration. Observation of STIRAP 
should be possible in devices within present fabrication 
standards, provided both design and operating condi- 
tions are carefully chosen. 

We have shown that efficient population transfer re- 
quires optimizations of the tradeoff between large enough 
pump coupling and the implied larger sensitivity to low- 
frequency noise. To this end the CPB should be biased 
slightly off-symmetry in a region where low-frequency 
fluctuations of the energy spectrum are linear in the fluc- 
tuations X of the control parameter qg . 

We have shown that the noise is conveniently analyzed 
by mapping it onto fictitious correlated fluctuations of 
the detunings (see Fig. [5]). This simple picture emerges 
because, despite of the complications brought by the mul- 
tidimensional space of parameters, the efficiency for STI- 
RAP is shown to depend essentially on noise channels rel- 
ative to the trapped subspace only. The relevant channels 
can be fully characterized by operating the nanodevice as 
a qubit, as in Refs. [29l - l33i 

We have found that the tradeoff is summarized by a 
single figure of merit, given in Eq.(|n|), which indicates 
favorable conditions for observation of STIRAP. Its re- 
markable dependence on features of the three-level spec- 
trum of the device (energy correlations, symmetries) sug- 
gests that band structure engineering may play a key role 
in determining optimal design solutions. This analysis, 
together with other already available tools, as improve- 
ments in materials and control circuits, besides a system- 
atic investigation of parameters and pulses crafting, guar- 
antees room for further improvement of the efficiency. 

In this work we did not consider other noise sources (as 
the readout circuit or critical current noise), which are 
possibly coupled to the device in channels "orthogonal" 
to the drive. This is because in the successful regime for 
STIRAP they lead to minor effects in CPB's^. They can 
be easily accounted for by a slight generalization of our 
approach, allowing for independent noise sources. Notice 
that each noise source could determine its own correla- 
tions of (5 and 8 p. 

We remark that the physical picture emerging from 
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this work applies to the whole class of superconduct- 
ing nanocircuits, used so far for implementing quantum 
bitSii^. Our full analysis applies to flux-qubits^ where a 
coordinate-parity selection rule holds^ and a symmetry 
point exists, except that two orthogonal noise sources 
(flux and critical current plus charge^) should be taken 
into account for accurate predictions. It also applies 
to phase-qubits^ where only linear fluctuations are im- 
portant^, but detunings are differently correlated. In 
all these devices the figure of merit analogous to that 
of Eq.® can be used to characterize the effect of low- 
frequency noise versus efficient coupling. 

A natural extension of our work is the investigation of 
dynamic diffusion for correlated phases/detunings in the 
experimentally relevant case of 1//" noise, and to which 
extent the dependence on the drive intensity of the re- 
silience to low-frequency noise can be used for some effec- 
tive dynamical decoupling. Moreover it has been pointed 
out that the presence of one or few more strongly cou- 
pled fluctuators may deteriorate the efficiency of ideal 
STIRAP— and it would be interesting to extend the in- 
vestigation to the LZ scenario. Finally circuit-QED2^ 
based architectures are natural candidates for the imple- 
mentation of STIRAP with quantum fields, the physics 
related to BBCN^ must be studied in this broader sce- 
nario. 



Appendix A: More on the CPB 

1. Driven three-level effective Hamiltonian 

Manipulation of the quantum state is performed by 
adding to the DC part of the gate voltage AC microwave 
pulses with small amplitude, qg — > qg + qc{t). The result- 
ing Hamiltonian can be written as 

Hit)=Hoiqg) + A{t)h (Al) 

where A{t) = —2Ecqc{t)- The effective three-level arti- 
ficial atom Hamiltonian is obtained by projecting H{t) 
onto the subspacc spanned by the three lowest energy 
eigenvectors 10^), i = 0, 1, 2 of i?o(<Zg) 

H{t) = Y,E,m{^^\+Ait) ^n,, 10,) (</),! (A2) 

i ij 

where — {(j)i\h\(l)j) . The STIRAP protocol can be 
carried out if we let A{t) = As{t) cosLOgt + Ap{t) cosujpt. 
We then perform the RWA, by retaining only quasi reso- 
nant off-diagonal and corotating terms of the drive, which 
simplifies to 

A(t)n -> HRwit) = I [ni2 A(i)e*"=*|0i)(02| ^^^^ 

+no2^p(t)e'"''*|(/)o)(02|] +h.c. 

Finally, the Hamiltonian is transformed to the doubly 
rotated frame, at angular frequencies ujg and Up via the 



transformation Urf = exp[i(w,j|</'i)(0i| + ^p\4>a){4'o\)t]- 
This yields an effective Hamiltonian H{qg) with the 
structure of Eq.(IT]), implementing the A configuration. 
Notice that — {(t)i\h\(j)j) play the same role of the 
dipole matrix elements in the definition Eq.([5]) of the 
Rabi frequencies. 

2. Charge-parity symmetry and selection rules 

Charge parity is a possible symmetry of wavefunctions 
in charge space which emerges because of the discrete na- 
ture of the momentum. Formally we introduce operators 
TLq = I? ^ which implement a refiection and 
then a translation in the charge space. If the parameter 
q is integer H^ always operates onto the same Hilbert 
state of discrete charges. It is easy to see that 

U-'[Y,in~qgf\n){n\]ll,^J2(''-q + qgf\n){n\ 

n n 

H-i[5]|n)(n±l|]H,=^|n)(nTl| • 

n n 

Therefore, one can seek for the invariance of the fam- 
ily of Hamiltonians ([3]). Symmetry points are found 
for qg — q/2, where HQ{qg) is invariant with respect 
to ^2qg- Since H^ = 1, for symmetric Hq eigenval- 
ues can be chosen with a well defined charge parity 
^2qg\4>j{qg)) = (—!)■' \4'j{'lg)) and parity selection rules 
hold such that for states of different parity charge matrix 
elements vanish, {(j)j\n\(f>i) = 0. 

Appendix B: Fluctuational behavior near the 
symmetry point 

We give a more detailed account on the effects of 
low-energy fluctuations close to qg = 1/2, displayed in 
Figs. lSllUl and on how they combine with high-frequency 
noise. On approaching qg = 1/2, fluctuations of the de- 
tunings turn from linear to quadratic. These fluctuations 
alone (thick red squares in the figures) would determine a 
nonmonotonic behavior of the efficiency on approaching 
the symmetry point. Indeed for 0.49 < < 0.495 fluc- 
tuations for X > yield smaller stray detuning 6{qg + x) 
than in the linear approximation, and the efficiency in- 
creases. However, approaching the symmetry point fluc- 
tuations (|(5|) exceed the linewidth 6i. Indeed, since this 
latter scales with flo (see sec. |V| and thus vanishes for 
qg 1/2, we find (|5|)/r2o oo and the efficiency should 
eventually vanish. However in this regime also the effect 
of ffuctuations of the couplings play a role. We study 
these fluctuations in linear and quadratic approximation, 
n-ijiqg + x) Ki nm{qg) + Aq2X + ^B^x^ ■ It turns out 
that only fluctuations of no2 are possibly relevant, and 
only in the regime where fluctuations of 5's are quadratic. 
However in this regime they spoil the picture based on 
fluctuations of (5's only. Indeed, for 0.49 qg ^ 0.495 



13 



smaller values of ^p{qg + x) for a; > compensate the 
positive effect of smaller S{qg + x). On the contrary, on 
approaching the symmetry point slow fluctuations of no2 
provide a nonvanishing coupling which is enough to yield 
a nonzero efficiency. Notice that this is true even at the 
nominal bias Qg = 1/2, where the selection rule is exact, 
since also in this limit {\S/no2\) is finite. 

In Figs. I8|10l we show the effect of linear (blue solid 
curve) and quadratic (blue dots) fluctuations of no2 . The 
two approximations differ for small J, indicating that the 
series expansion is likely not accurate enough. This is not 
a problem for our description of STIRAP, which in this 
regime is anyway suppressed by spontaneous decay. For 



larger values (J = 1.32, 2) the series expansion is seen to 
be accurate. 
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